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1. INTRODUCTION 


This report presents the results of a feasibility study of using magnetic 
torquing for Control Moment Gyro (CMG) desaturation and momentum dumping 
for the Apollo Applications (AAP) mission. In magnetic torquing, a torque 
is developed by placing a magnetic dipole moment in a magnetic field. This 
torque is applied to the vehicle, and its effect is counteracted by the CMG's. 
The result is to change the components of momentum in the Control Moment 
Gyros. The effectiveness of this technique to meet the requirements for 
CMG momentum dumping is the problem to be analyzed. The necessary control 
laws for such a system are derived, and a simple control system design is 
presented. The basic equations for designing either coreless or ferrite core 
magnetic torquers and representative design calculations are presented in 
a "tutorial format" to aid in understanding the problem. 

This study investigates an alternative to the gravity gradient technique 
for CMG desaturation because the gravity gradient technique has some undesir- 
able features and limitations. Since gravity gradient methods are not a 
part of this study, these problems will not be discussed and no comparison 
will be made between these two systems. 

A review of the Marshall Space Flight Center (MSFC) and Belcomm reports 
on the magnetic torquing performed as a part of the Controls Analysis sub- 
task 7 on AAP task KM205 is included in this report. 
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2. CONTROL LAWS 


2.1 VECTORING OF THE MAGNETIC DIPOLE 

To implement p. magnetic torquing system, control laws must be derived 
to vector a magnetic moment. In this section the physical equations for the 
control laws will be explained and a control system adequate for CMG momen- 
tum dumping proposed. 

A magnetic torquer, whether coreless or ferrite, may be represented by 
a simple magnetic dipole moment for the purposes of illustration and calcu- 
lation (References 1, 2). The basis for the system is the phenomenon exhib- 
ited by a magnetic moment in the presence of a magnetic field. A torque 
vector, T, will be exerted on a magnetic moment vector, m, i - a magnetic 
field vector, B, according to the following relation: 

T = m x B (1) 

The momentum to be dumped is represented by the vector H with compon- 
ents 1^ and H c (Figure 1). By the definition of the cross product, Equation 
(1) indicates that a torque cannot be generated collinear with the earth’s 
magnetic field. Consequently, the vector H c represents momentum which cannot 
be dumped using torque derived from the earth’s magnetic field vector with 
the instanteous value shown in Figure 1. Therefore, T d represents the torque 
required to dump the momentum H n in some given time interval. This propor- 
tionality yields the relationship 

T. = -KH (2) 

d n 

where K is a time dependent proportionality constant yet to be determined. 
Given a constant required torque, T, the proportionality constant is simply 
the reciprocal of time. 

The value of the total dump time t is specified by the total dark side 
flight time and some estimate of favorable B field orientation during this 
time. It will be obvious in a later section that the optimum power expendi- 
ture is effected if a minimum required m is produced. To accomplish this 
the constant K is set equal to 1/t. 
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B earth's magnetic field 

H momentum to be dumped 

M magnetic moment which may be 
generated in any direction 

torque desired 

H momentum colinear with earth's field 
c 

momentum perpendicular with earth's field 
Figure 1. Torque Vector Diagram 


This system may take the form of an iterative process as follows: 
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(4) 


where H x , Hy, and H 2 are the total momentum to be dumped. These values are 
determined as an average CMG momentum buildup during one orbit. The AH , 

X 

AHy, and AH^ terms represent instantaneous momentum to be dumped; T^, T^, 
and T z represent the desired torque to be produced. The terms t Q and t^ are 
the beginning and end times for torquing while t represents the instantan- 
eous time. 
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Since H^» and are continually being incremented, they are 

cime dependent. In order to simplify the calculations, the torque components 
V T y* and T z should remain constant. As t approaches t f , AH should ap- 
proach zero to allow a finite solution for the torque given by Equation 14. 
However, errors in calculations for AH may require limiting the term (t f -t) 
to be a positive nonzero quantity. At several iterations prior to the 
termination time, this limitation will be imposed. 


From Figure (1) it can be seen that 


and 


H„ 53 (H » B) B 

c ~ 9 

B- 

H = H - If 

n c 


(5) 

( 6 ) 


Substituting equations (5) and (6) into Equation (2) and using a vector 
identity 


T , = “ (B 2 H - (H • B) B) 

B 

T d = ~ (H X B) X B 


(7) 


Relating this equation to Equation (1) and solving for m 

in X B = ~ (H X B) X B 

B 

m - 4- (H X B) (8) 

B" 

From Equation "3) , the optimum magnetic nrnent vector dumps the mo- 
mentum H, given the earth’s magnetic field B. The instantaneous representa- 
tion of Figure (1) is somewhat misleading since the orientation of the 
earth's magnetic field with respect to the momentum to be dumped will change 
along the orbital path, the momentum H c will no longer be colinear with the 
earth’s field. The capability to dump this momentum will then be acquired. 

There are two alternatives to the problem of proper magnetic moment 
vectoring. The first method involves the use of three orthogonal torquing 
units. By creating a magnetic moment of the correct magnitude in each 
torquer, the vector sum can be adjusted to produce the desired moment vector. 
In the second method a single torquer is gimballed to the proper orientation 
which creates a moment of the necessary magnitude. The gimbal structure 
must have two degrees of freedom so that the magnetic dipole can be vectored 
in any arbitrary direction. 
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The development of a control law for a system of three orthogonal tor- 
quers is straightforward. The desired magnetic moment m is determined by 

the use of Equation (8). The magnitude of the magnetic moment produced by 
the torquer will be a function of the applied current !. In a coreless 

torquer this relationship is linear; while in ferrite torquers, it may be 
nonlinear. The moment may then be represented by 

m « Cl (9) 

where C Is a proportionality constant determined by the physical character- 
istics of the coil and 1 is a vector representing the vector sum of the 
three orthogonal moment producing currents. In the case of a coreless tor 
quer, for instance, this would equal the turn-area product: 


where 


C = NA 

N * number of turns in coil 
A •= area of coil 


Combining Equations (8) and (9) and expressing the cross product in scalar 
form, the desired relationship is expressed in Equation (10). 
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The AH terms represent the instantaneous momentum to be dumped. Since 
B, AH and K will be functions of time, toe current in each torquer will also 
be a function of time. Except for the deterministic constant K/C, the cur- 
rent inputs to the torquers are specified. 


If a single gimballed torquer method is used, an equation of the same 
form as Equation (10) results. In this case however, the magnitude of the 
magnetic moment is required. Thus 
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This equation is again time dependent. From this, the necessary relation- 
ships follow: 

1 m [ - (m ^ + m ^ (12a) 

1 x y z ' 

* jmj (12b) 

c 

The constant C again depends on the specific torquer used, and I is the 
current requirement. 

Figure (2) represents a single torquer in a two gimbal system. Two 
degrees of freedom allow the magnetization vector to be properly 



positioned. Transformation from the two gimbal system shown in figure (2) 

to body coordinates is accomplished by the use of Equation (13) where 0^ and 

I, are the outer and inner gimbal angles respectively. The transformation 
8 

takes the form: 
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it therefore follows that 


From Equation 
components of 
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(13) the gimbal angles may be determined from the specified 
required magnetic moments in body coordinates: 
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2,2 DETERMINATION OF THE EARTH’S MAGNETIC FIELD 

The earth's magnetic field vector may be determined in several ways. 

The geomagnetic field may be stored in the onboard computer in the form of 
a table, with the independent variables identified as longitude and latitude 
and the dependent variables as the components of the earth's field in a 
local vertical reference frame. A fixed altitude would be assumed. This 
method would require a two-dimensional (latitude - longitude) interpolation. 
The stored table may require a storage of approximately 300 fixed values 
for the earth’s magnetic field and an undetermined amount of storage ca- 
pacity for the interpolation routine. Another possible method is to use a 
spherical harmonic model for the geomagnetic field. Again large storage 
capacity may be required to store the fundamental harmonic equations and 
coefficients. As a third alternative, three orthogonal magnetometers could 
be used as sensors to provide the scalar components of the field vector in 



body coordinates. The software requirements for this method are much sim- 
pler than the other two mentioned, requiring only that the magnetometer 
measurements be inputs to the onboard computer. However, the major problem 
with this method is to provide adequate isolation of the additional hard- 
ware, namely the three magnetometers, from the magnetic torquers. 

Any one of these methods appears feasible as a method of determining 
the earth’s field using onboard computations or sensing. However, the 
tradeoffs among these methods is beyond the scope of this investigation. 
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3. GENERATION OF MAGNETIC MOMENTS 


3.1 ASSUMPTIONS 


Several assumptions are necessary before the design criteria are dis- 
cussed so that a complete understanding of the design is achieved. To limit 
the complexity of the calculations* the value of the earth's magnetic field 
is taken to be constant at a value of 0.3 gauss. This value represents the 
average magnetic field magnitude at 230 nautical miles based on a dipole 
model of the earth's magnetic field, A momentum value of 1000 foot-pound- 
seconds is taken for the CMG dump requirement and it is assumed that one 
half orbit (45 minutes) is available for dumping. This dictates the torque 
requirement of 0.5 Newton meters necessary to dump accumulated momentum. 
These values are assumed constant throughout. 


3.2 CORELESS TORQUERS 


The coreless torquer is considered as a current loop. Conceptually, 
the force acting on a current loop is caused by the magnetic dipole moment 
of the loop, but, in actuality, it results from the force acting on a cur- 
rent element: 


where 


and 


F = q (V X B) = I X B 

F represents force 
q represents charge 
v is the velocity of the charge 
i represents current 


(15) 


The torque acting on the loop may be developed from this equation. Figure 
(3) is a representation of a current loop. Referring to Figure (3), Equa- 
tion (15) may be written in the form: 


dF = Id l X B 


(16) 
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Figure 3. Current Loop 

where I is the current flowing in the loop and dJl is an incremental element 
of the loop. From the definition of torque, 

d T = r X df (17) 

From figure (2) the magnitude of 

d£ = Rd6 (18) 

and r = Rsin 6 (19) 

Combining Equations (16), (17), (18), and (19), the resulting magnitude 
equation represents the present system: 

dT = IR 2 B Sin 2 0 d 6 (20) 

Integrating this equation and evaluating the integral, the resultant torque 
is 

T = 7T R 2 I B (21) 

Substituting A = tt R 2 into the equation and I = Ni, a more familiar form 
results : 

T = A Ni B (22) 

where the product ANi is defined as the magnitude of the magnetic moment m. 

Based upon the assumed values for T and B, Equation (22) may now be 
used to calculate the magnetic dipole moment necessary from the coreless 
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torquer , and a value for m of 0.166 X 10 amp - m is determined. Since 
the magnitude of m is equal tc the product ANi, it is obvious that some 
tradeoffs are available in the design of the loop. Since A * rR , the mag- 
netic dipole moment is a function of the radius R squared, and it is appar- 
ent that the larger the loop’s radius, (and thus the are?) the lower the 
current requirement. The power requirement is solely a function of the 
ohmic loss and, therefore, proportional to the radius of the loop and the 
square of the current in the loop. Larger loops, consequently, dissipate 
less power than smaller loops. By means of the same reasoning, it may be 
shown that smaller loops have more mass for fixed torque requirement. The 
various tradeoffs between radius, power, and weight are indicated in the 
curves given in Figure 4. 



LEGEND 


CURVE 

RADIUS 

(fO 

NUMBER 
OF TURNS 

WEIGHT 

(lb) 

A 

1.6 

2 X 10 4 

67 

B 

3.3 

5 X 10 3 

33 

C 

10.0 

6 X 10 2 

25 

D 

1.6 

2 X !0 6 

670 

E 

3.3 

5 X 10 5 

334 

F 

10.0 

6 X 10 3 

246 


ASSUMPTIONS 

• B = 300 MILUGAUSS (CONSTANT EARTH MAGNETIC FIELD) 

• DUMP TIME = 45 MIN 

• WINDINGS CONSIST OF NUMBER 30 WIRE THROUGHOUT 


Figure 4. Coreless Magnetic Torquer Performance 
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Example calculations will be performed for a 2 meter and a 10 meter 
diameter loop. The windings consist of number 24 annealed copper wire with 
an electrical resistance of 84.21 ohms per kilometer and a mass of 1.82 
kilograms per kilometer. A current value of 0.10 ampere will be used in 
each case. 


In the case of the 2-meter diameter loop, 

NAI = 0.166 X 10 5 amp - m 2 

A * 7tR 2 * 3.1416 rn 2 

Therefore, the number of turns necessary is 

N * 5.29 X 10^ turns 


The total length of wire, 4, follows: 

i = N2ttR * 33.25 X 10^ meters 


The resistance of the coil is 28.0 kilohms and therefore the power 
dissipated is 


P dis 

The weight of the loop is 


280 watts 


W = 605 kilograms = 1320 pounds 

In the case of the 10 meter diameter loop, the area is given by 

A - ttR 2 = 78.5 m 2 

The number of turns necessary in the coil is 

N = 2115 turns 


The total length of wire necessary, £ , is 

3 

£ = N2ttR = 66.4 X 10 meters 


The resistance of the wire comprising the coil is 5,590 ohms, and the 
power dissipated in ohmic losses is 


P dis 

The weight of the loop is 


55.9 watts 


W = 121 kilograms (266 pounds) 
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These calculations indicate the obvious tradeoffs which can be made 
between the weight, radius, and power parameters. For example, if weight 
is of more concern than power consumption, lower weight can be achieved 
through smaller diameter wire. For the 10-meter air loop, a number 30 
annealed copper wire will decrease the weight to 30.2 kilograms (66.5 pounds) 
with a power consumption of 225 watts. 


Although the derivation will not be carried out here, it may be shown 
(Reference 2) that the magnitude of the magnetic field resulting from a 
current loop, along the axis of the loop is 


B 


1 

2 


(Z^ + a z ) 


2x2/3 


(23) 


where a is the radius of the loop and z is the distance from the center of 
the loop along the axis. Substituting the appropriate values, the magnitude 
of the magnetic field in the plane of the loop, for the 10 meter diameter 
loop will be 

B * 0.534 gauss 

At a distance of one meter, this value is 

B * 0.502 gauss 

The field strength is only about twice that of the earth’s field and rep- 
resents no difficulty to any electronic equipment onboard. 

3 . 3 FERROMAGNETIC MATERIAL 

All materials exhibit some magnetic properties. Depending on their 
magnetic behavior, substances can be classified as diamagnetic, paramagnetic, 
and ferromagnetic. If a magnetic field is applied to these materials, the 
magnetization in the diamagnetic materials is opposed to the applied field. 

In the paramagnetic materials, the magnetization is in the same direction 
as in the applied field. In both cases, however, the effect is weak. Ma- 
terials in the ferromagnetic group, on the other hand, exhibit very strong 
magnetic effects. It is these materials which will be considered for 
addition as a core to the torquer with the expectation of obtaining some ben- 
efit from their inclusion. 

The characteristic of ferromagnetic material is to produce a magnetiza- 
tion vector proportional to the field applied to it. This magnetization 
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vector may be viewed as a volume aggregate of atomic magnets and is defined 
as the magnetic dipole moment per unit volume. It is then possible to de- 
fine the magnetic dipole moment of a volume of ferromagnetic material as 
the product of the volume of the material v and the magnetization vector 
M. Thus, 

m * Mv (24) 

For a given magnetization vector and specified magnetic moment, the required 
volume of material can be determined. The magnetization vector will be a 
function of the material considered and the magnitude of the applied field. 

Two distinct groups of ferromagnetic material are considered for use 
as torquer cores. The first group (high u ferrites) is represented by two 
materials, Supermalloy and Vanadium Permendur. In general, these materials 
exhibit a high permeability with little magnetic retention; that is, when 
the applied field is reduced to zero, the magnetic field within the material 
returns to approximately zero. The second group (the permanent magnetic 
ferrites) is represented by Alnico V. This material exhibits a high per- 
meability as well as the ability to retain a strong magnetic field within 
the material as the applied field is returned to zero. This ability is 
called magnetic remanence, and its value is dependent upon the material used 
and the past history (hysteresis) of applied magnetic fields. Pertinent 
data for these materials are listed in Table 1. 


Table 1. Ferromagnetic Material * 


Name 


Br 

remanence , 
(gauss) 

Maximum perme- 
ability (gauss/ 
oersted) 

He Coer- 
cive Force, 
(oersted) 

Bs=AttMs Satur- 
ation Indi- 
cation (gauss) 

Alnico V 


12,500 

15 

660 

nil 

Supermalloy 


nil 

1,000,000 

0.002 

7,900 

Vanadium Permendur 

nil 

4,500 

2.0 

24,000 


* 

See Reference 

5 




3.4 DEMAGNETIZATION 

With the introduction of ferromagnetic material, it is necessary to 
discuss the concept of a demagnetizing field. This demagnetizing field 
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accounts for a noncontinuous ferromagnetic path (Reference 5). When a rod 
is magnetized by an applied field H , the ends carry magnetic poles which 
cause magnetic fields in all parts of the rod. Normally, these fields are 
directed in the opposite direction to the applied fields and are, therefore, 
called demagnetizing fields. The true field acting on a given section of 
the bar,(e.g., its middle) is then the resultant of the applied field and 
the demagnetizing field AH 

H = H - AH (25) 

a 

The demagnetizing field is approximately proportional to the magnetiza- 
tion 

AH = NM (26) 

where N is the demagnetization factor. It is apparent from Equation (26) 
that a large demagnetization factor will give a large demagnetizing field. 
For the ferrites, the demagnetization factor is a function of the geometry 
considered. Equations to determinn demagnetization factors have been pub- 
lished. Tables and charts based on the solutions of these equations for 
cylinders and ellipsoids of revolution have also been published (Reference 
5, 6, 7). Table 2, for example, lists demagnetization factors for various 
length to diameter (J./D) ratio cylinders. 

Table 2. Demagnetization Factors for Cylinders 


L/D Demagnetization Factor 


1 

0.27 

2 

0.14 

5 

0.04 

10 

0.0172 

20 

0.00617 

50 

0.00129 


A large L/D ratio for the solenoid core and a corresponding low value 
of demagnetization factor are desirable. How large L/D may be is a function 
of space limitations. For the purpose of calculations, a ratio of 20 is 
chosen and the resulting demagnetization factor is 0.006. 
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3.5 SOLENOIDAL TORQUERS 

A helical coil, or solenoid, is often used to produce a magnetic field. 
The coil would consist of N turns of wire carrying a current I. Although 
the proof will not be carried out here, it may be shown (Reference 4) that 
a torque is exerted on a solenoid in an external field B so that 

T * NIAB (27) 


where A is the area of a cross section of the solenoid. This is the same 
expression derived for the loop (Equation 22, Section 3.2). Torque produced 
for the loop is a limiting case of that for the solenoid. A ferrite core 
was not considered for use with the loop because the loop is an inefficient 
method to produce a field in a long cylinder due to the field pattern of 
the loop. If a thin cylinder is considered, restricting it to the thickness 
of the loop and its location to the center of the loop, the demagnetization 
factor is approximately unity (1), and, therefore the ferrite is ineffective. 

The magnetomotive force at the center of a long solenoid i.e., where the 
length of the helix is much larger than the mean radius of the helix, can be 
shown (Reference 4) to be 

H = NX (28) 

a 

where a is the radius of the loop. The field necessary to produce the de- 
sired magnetization vector is given by 


H * H_ + AH 

a 

where H represents the field as seen b> the ferrite, AH represents the 

ci 

demagnetizing field and H is the fxeld applied to the ferrite. The magnetic 
field is given by the relationship 


B, 


U(H)H 


ferrite 

where y(H) is the nonlinear relationship known as permeability. 


The solenoid is used to apply a magnetic field to a cylinder of fer- 
rite material. It will be wound on the surface of the cylinder and be re- 
stricted to this area. Calculations of the physical parameters of the 
torquer will be accomplished for each ferromagnetic material. As was dis- 
cussed earlier, a torquing system using the solenoidal torquers may consist 
of three stationary orthogonal torquers or one gimballed torquer. 
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3.5.1 Supermalloy 

Consulting the magnetization curves for Supermalloy (Reference 5), 
an operating point must be chosen. In this instance, 

M - 557 X 10 3 amps/m 

was chosen below saturation to avoid Inefficient utilization of the ferrite 

Therefore, with the desired torque and assumed magnitude of the geo- 
magnetic field, the volume of material necessary is 

V - 0.0299 m 3 

With an L/D ratio of 20, this results in the radius of the cylinder equal 
to 0.619 meter (2.43 inches) and the length of the cylinder equal to 2.479 
meters (97.4 inches). The resultant demagnetizing field using a demagnetiz 
ing factor of 0.006 is 

AH = 3.34 X 10 3 amps/m 

From the magnetization curves, a resultant field of 0.06 oersted or 4.78 
amps/m is necessary to produce M. Therefore, the total applied field 
necessary is the sum 

H = H + AH * 3.344 X 10 3 amp turns 

a 

The physical quantities of a solenoid are related to the field by Equation 
(28). Now, since the solenoid is woun$ uniformly over the cylinder, and 
the cylinder is 2.479 meters long, 

NIL * nl - 8.28 X 10 3 amp turns 

where n equals the total number of turns. If the current I is 0.5 ampere, 

then ^ 

n * 16.56 X 10 turns 

and the length of wire necessary for a cylinder of this size is 

Z = 2-nrn = 6.44 km 

Using number 24 annealed copper wire, the coil will weigh 11.75 kg (25.9 
pounds) and have a resistance of 544 ohms. The power loss will then be an 
I^R loss and amounts to 136 watts. The total weight of a single torquer 
will consist of the core weight of 573 pounds and the coil weight of 25.9 
pounds, totaling 598.9 pounds. There is virtually no hysteresis loss en- 
c ‘"erec with Supermalloy. 
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3.5.2 Vanadium Permendur 


For Vanadium Permendur, a point is chosen on the magnetization curve 

where , 

M ■ 1.595 X 10 amps/m 

This point, far below saturation, was chosen to take maximum advantage of 
the slope of the magnetization curve. The volume of material required is 

3 

v ■ 0.0104 m 

and, for an L/D ratio of 20, results in a cylinder with a radius equal to 
0.0436 meter (1.74 inches) and a length of 1.746 meters (68.5 inches). 

The resultant demagnetizing field using a demagnetization factor of 0.006 
is 

AH * 9.55 X 10^ amps/m 

From the magnetization curve a resultant field of 

H * 1.2 X 10^ amps/m 

d 

is necessary, and, therefore, the applied field consists of the sum 

H * + AH » 10.75 X 10^ amp-turns 

The cylinder, being 1.746 meters long, requires an amp-turns figure of 

NIL * nl * 18.75 X 10^ amp-turns 

For a current value of 1 ampere 

3 

n ■ 18.75 X 10 turns 

and a length of 5.11 kilometers of wire is necessary. Again using number 
24 annealed copper wire, the weight of the wire is 9.37 kilograms (20.65 
pounds) and the resistance 434 ohms. Power dissipation in the wire is 434 

3 

watts. Vanadium Permendur has a discernable hysteresis loss of 600 joules/m 
and thus, an additional power dissipation of 6.26 joules (watt/sec) is 
suffered during each turn-on and turn-off sequence. The weight of single 
torquer, solenoid windings and core material, amounts to 94.97 kilograms 
(209.15 pounds). 

For the purpose of comparison, if the L/D ratio is now changed to 50 

r 

for a Vanadium Permendur torquer, the cylinder would have a radius of 1.46 
inches and a length of 12.4 feet. Because of the increased L/D ratio, the 
demagnetization factor decreases by a factor of 4.8 to 0.00129. Recalculat- 
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ing the pertinent quantities, if the current is decreased to 0,5 ampere, the 
power dissipation in the wire is now 106.5 watts with the same torquer 
weight of 209.15 pounds. This is a substantial savings in pever. 

3.5.3 Alnico V 

The permanent magnetic material, Alnico V, exhibits the phenomenon known 
as remanence. When a field of sufficient magnitude to b^ing the material into 
or near saturation is applied and subsequently removed, a strong magnetization 
vector will be retained within the material. Saturation in this material is 
represented by a inductance of 

B * 12,500 gauss 

Repeating the type of calculations carried out above 

M * 962 X 10 3 amps/m 

AH * 5.772 X 10^ amps/m 

Therefore, the required applied field is 

* H = AH + = 77.46 X 10 3 amps/m 

The volume of required material is 

V = 0.01722 m 3 

and therefore, for an L/D ratio of 20, the cylinder has the dimension of 

r = 0.0515 meters (2.025 inches) 

L = 2.102 meters (83.6 inches) 

For these values, the required coil will consist of 

3 

nl = 162.5 X 10 amp-turns 

Since this is an exceedingly large number, a current of 5 amperes is chosen 
to avoid astronomical coil weights. A coil wire length of 10.52 kilometers 
is necessary with a resistance of 886 ohms and a weight of 19.15 kilograms 
(42. 2 pounds). The power dissipated through ohmic losses is 22.2 kilowatts. 
The hysteresis loss for this material during a torquing on-off cycle is 
approximately 330 joules. The total torquer weight is 144.75 kilograms 
(319.2 pounds) . 

V 

4 Although the power requirement is high, this represents the instantan- 

eous peak power necessary. If a current pulse of finite width is used, as 


3-11 



is possible with a permanent magnetic material, the average power over a 
single torquing period will be much lower. If, for example, a 6-second 
turn-on and turn-off pulse is used, an average power requirement of 50 
watts is necessary. 

3.5.4 Tradeoff Factors 

The amount of material necessary for a torquer is a constant for a 
specific material. The weight and power requirements of the solenoidal 
coil however, are variable, and tradeoffs are possible. For example, if a 
number 20 annealed copper wire is used with the Vanadium Permendur core 
instead of the number 24 wire, with a L/D ratio of 50 at a current of 1.0 
ampere, the power dissipation will be 85.5 watts, and the coil will weigh 
26.2 pounds. Also, for a given wire size the power-weight product is con- 
stant; for example, if the current is decreased by a factor of 10, the 
power will decrease and the weight will increase, each by a factor of 10. 

3.5.5 Field Quantities 

The magnetic field exterior to the torquer is an important quantity 
because of instrumentation onboard the vehicle. For a sufficiently long 
solenoid, the field along the side of the solenoid will be negligible 
(Reference 2). The only field of consequence will be from the end faces. 



Figure 5. Solenoidal Torquers 

In Figure (5), the expression for the magnetic field at point P along the 
axis of the solenoid will be 

B = ^ NI (cos 02 ~ cos 0^) 

B = y p M 
2 


& + x 


\/r 2 +(£J-x) 2 


«_1 
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Carrying through calculations for Vanadium Permendur with 

x * 1.0 meter 

r = 0.0436 meter 

i = 1.746 meters 

2 

Bf = U Q (H 4- H) * 2.02 webers/m 

we have : 


B = B i 
2 

then the field at one meter is 


SL + x 


y r^+ (£+x)^ 


+ x 


B (1 meter) = 1.01 


2.746 


/ 0.0019 + 7.5405 v/0.0019 + 1 

- 8.484 gauss 


Figure 5 is representative of the solenoidal torquers with ferrite cores as 
a group. At three meters, the same calculations show that the field is re- 
duced to 0.404 gauss or slightly more than the magnitude of the earth's 
field . 


3.6 SPHERICAL TORQUERS 

A spherical core wound with three orthogonal windings generates a general 
magnetic dipole with the advantage of a weight and space savings over sys- 
tems using three separate cores. Although the weight and space saving seems 
attractive, the demagnetization factor of 0.333 for the sphere represents 
a distinct disadvantage. The implications of this will become more obvious 
after later calculations. 

The spherical torquer will be excited with a coil wound on the sphere 
with uniform axial density (Reference 1). Referring to Figure 6, this will 
be represented by 

K = i — - m - sin 6 (30) 

4ttR J 

The field inside the sphere will have the form 

H - — — (i r cos 6 t L sin 0 ) (31) 

1 2TTR 3 
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and outside the sphere, the field will be that of an ideal magnetic dipole 


H = ~nr~ (I cos 6 + i fi h sin 0 ) (32) 

2ttR J r U 



Figure 6. Spherical Torquer 

Using Supermalloy, due to its extremely high permeability and negligible 
remanence field, 

3 

M = 557 X 10 amps/m 

Now, with a demagnetization factor of 0.333, the demagnetizing field is 
equal to 

AH = 185.481 X 10^ amps/m 

With a necessary field of 4.78 amps/m, the applied field must now obtain the 
value 

H = 185.486 X 10^ amps/m 

The coil windings add a significant torque as a result of the magnetic di- 
pole moment they produce. This torque results from the large demagnetiza- 
tion factor (0.333 for the sphere versus 0.006 for a L/D = 20 solenoid). To 
calculate the dimensions of the material necessary, it will be noted that 
the total magnetic moment from the core and the winding takes the following 
form 

Mv + NIA = m (33) 


3^14 



where v is the volume and A is the area of the sphere. The coil necessary 
to produce the necessary magnetic field is 


NI = 3m_ = H • 3R 

2 

2vR 


or approximately 


NI = M • R 

Equation (33) then takes the form 

= 3m 
16 ttM 

Substituting the appropriate values and carrying out the calculations the 
value for the radius is obtained: 


R = 0.1694 meters = 4.48 inches 

The volume of the sphere is then 

V = 0.02036 m 3 

This corresponds to a weight of 178.51 kilograms (394 pounds) for the core 
material. 

The coil necessary is now 

NI = H * 3R 

= 3 (185.486 X 10 3 ) (0.1694) 

= 94.29 X 10 3 amp-turns 

Since the same tradeoffs among weight and power made earlier for the solen- 
oid are also made here, a current equal to 1.0 ampere is chosen. To calcu- 
late the necessary length of wire for a sphere, the average radius must be 
known. This is determined as follows: 


R 

av 


— f Rm 2 
Rm o m 


- x ) dx 


Rm 


x 

2 


/■Rm 2 - 


Rm 


2 2-1 
x + sin , x v 

2 '‘Rm } 


= Rm (0.7854) = 1330 meters 


3r<L5 




Figure 7. Radius Calculation 


«r 

% 

4 


Then 


i = 2ttR N = 78.82 kilometers 

av 


Again using number 24 wire, the resistance of the wire is 6.448 kilohms 
and the weight of a single coil is 143.3 kilograms (316 pounds). The 
power dissipation is now 

P , . = I 2 R = 6448 kilowatts 

dis 

The total weight of the torquer for one axis momentum generation capability 
is 321.8 kilograms (710 pounds). 

To check the accuracy of these calculations, the magnetic dipole moment 
must be calculated for the core and the windings. 

m = Mv 

C 3 

= (557 X 10 J ) (0.02036) 

= 0.05259 X 10 5 amp - m 2 

Summing these: 

m + m = 0.16613 X 10“* amp - m 2 

c w 

= m as desired 


where m^ equals the magnetic moment of the core and m^ equals the magnetic 
moment of the windings. It is fairly obvious by looking at the physical 
parameters for the windings that construction of this unit is virtually im- 
possible. If the amount of wire necessary is reduced by a factor of 10 
to a more reasonable figure, the power requirement increases to 64.48 kilo- 


watts . 
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As can be seen from Equation (32), the form of the magnetic field is 
that of a magnetic dipole. Taking the derivative with respect to 0 of 
Equation (32) and setting the result equal to zero, the maximum value of the 
external field occurs at an angle of approximately 26 degrees from the mo- 
ment axis. At this angle and using the above values, the magnetic field 
at a distance of one meter from the surface of the sphere was calculated to 
be 

B (1 meter) = 23.4 gauss 

This value is a factor of 2 larger than the field values for the solenoidal 
torquers 

3.7 TRADEOFF DISCUSSION 

Three orthogonal coreless loops wound around the S-IVB represents the 
minimum weight and power configuration. The magnetic field vector developed 
by the torquer has a magnitude approximately twice that of the earth's mag- 
netic field at the center of the loop. This appears to be at a safe level 
although no data are presently available on the magnetic susceptability of 
onboard equipment. This is one area which must be further investigated 
before final qualifications can be considered. If the coreless loops are 
not acceptable because of space, the alternative methods are represented by 
the ferrite torquers. The magnetic field produced by the ferrite torquers 
is many times the magnitude of the earth’s magnetic field, and its accept- 
ability is questionable. Considering only weight and power, a tradeoff com- 
parison can be made using the results presented in Table 3. 
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4. CONCLUSIONS 


Other reports issued by Marshall Space Flight Center , Bellcomm and NASA 
have documented the results of analyses made to determine the possible uses 
of magnetic torquing. Most of these reports are concerned with continuous 
attitude stabilization for AAP . However, Marshall Space Flight Center pub- 
lished a memorandum dated November 8, 1968, in which CMG desaturation is 
considered (Reference 9). This report concludes that magnetic torquing 
for CMG momentum dumping should be considered. The Bellcomm letter dated 
February 11, 1969 (Reference 10) states that the feasibility of the concept 
of attitude stabilization has been established. The Marshall report which 
concludes that the magnetic torquing method is feasible, is the only analysis 
dealing with the use of magnetic torquing exclusively for CMG desaturation. 
Discussions presented earlier in this report support this conclusion within 
certain stated qualifications. However, to arrive at an unqualified conclusion 

a more detailed investiga- ion in the form of a simulation is necessary to 
determine dumping capability. It appears that Bellcomm is pursuing this, 
although no results have as yet been published. Weight and power figures 
for magnetic torquers are included in Table 4 as published by Marshall and 
Bellcomm. They are included here for the purpose of comparison. Since no 
indication of computational assumptions accompanied the Marshall and Bell- 
comm results, it is not possible to do more than note the discrepancies 
between power and weight when compared to those given Table 4. 


4-1 



TABLE 3: COMPARISON BETWEEN CORELESS AND FERRITE CORE TORQUERS 
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Earth’s Magnetic Field Equal To 3.0 x 10 Webers/M 


TABLE 4. COMPARISON DATA FROM OTHER ANALYSES 


MARSHALL SPACE FLIGHT CENTER REPORT (REFERENCE 9) 


CORF. 

TYPE 

COIL 

WEIGHT 

CORE 

WEIGHT 

TOTAL 

WEIGHT 

CONTROL'* 

TORQUE 

POWER 

Aire ore 
3 . 5 b m d i a . 

210 lbs 


210 lbs 

. 5 Nm 

400 watts for 30 min 
1000 Nms 

Same 

53 .o 


53.6 lbs 

.16b Nm 

1 50 watts for 40 min 
1000 Nms 

A i rcore 
h . b m d i a . 

17.4 lbs 
92 . Z lbs 


17.4 lbs 
52.3 lbs 

.166 Nm 
.166 Nm 

150 watts for 40 min 
1 000 Nms 

50 watts for 90 min 
1000 Nms 

A i n i c o ^ 
4.5 in d j a 
90 in long 

1 lb 

356 lbs 

3 57 lbs 

. 5 Nm 

16 watts for 30 min 
1000 Nms 

A 1 n i c o 5 
3 i n d i a 
bO in 1 ong 

Torque de 

1 lb 

ve loped in 

116.5 lbs 
a magnetic field 

117.5 lbs _ 
of 3.7 X 10" ’ 

.166 Nm 

? 

W/M“ 

2.42 watts for 90 min 
1000 Nms 


BELLCOMM (Reference 11) 


CORE 

TYPE 

COIL 

WEIGHT 

CORE 

WEIGHT 

TOTAL 

WEIGHT 

CONTROL 

TORQUE 

POWER 


15 ft long 
12.8 in d i a 

unknown 

unknown 

7860 lbs 

unknown 

Peak power 52C 
Average power 

) watts 
3 1 0 watts 

8 f t 1 ong 
3-1/8 in dia 

unknown 

unknown 

233 lbs 

unknown 

Peak power 48 
Average power 

watts 
31 watts 

5 ft long 
2-3/8 in dia 

unknown 

unknown 

91 lbs 

unknown 

Peak power 31 
Average power 

watts 
20 watts 
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